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ABSTRACT 

Because the opacity of clouds in substellar mass object (SMO) atmospheres depends on the 
composition and distribution of particle sizes within the cloud, a credible cloud model is essential 
for accurately modeling SMO spectra and colors. We present a one-dimensional model of cloud 
particle formation and subsequent growth based on a consideration of basic cloud microphysics. 
We apply this microphysical cloud model to a set of synthetic brown dwarf atmospheres spanning 
a broad range of surface gravities and effective temperatures (g SU rf = 1-78 x 10 3 - 3 x 10 5 cm 
s~ 2 and T c ff = 600 - 1600 K) to obtain plausible particle sizes for several abundant species (Fe, 
Mg 2 Si04, and Ca2Al 2 Si07). At the base of the clouds, where the particles are largest, the particle 
sizes thus computed range from ~5 /jm to over 300 ^m in radius over the full range of atmospheric 
conditions considered. We show that average particle sizes decrease significantly with increasing 
brown dwarf surface gravity. We also find that brown dwarfs with higher effective temperatures 
have characteristically larger cloud particles than those with lower effective temperatures. We 
therefore conclude that it is unrealistic when modeling SMO spectra to apply a single particle 
size distribution to the entire class of objects. 

Subject headings: atmospheres: clouds, condensation, grains: fundamental parameters — stars: low 
mass, brown dwarfs, substellar mass objects, L dwarfs, T dwarfs, spectroscopy, atmospheres, spectral 
synthesis 



1. Introduction 

Substellar mass objects (SMOs) arc fundamen- 
tally more complex than Sun-like stars because 
of the formation of molecules in their cool outer 
layers. These molecular species cause the spec- 
tra to deviate strongly from blackbody values in a 
multitude of spectral bands (Leggett et al. f999; 
Kirkpatrick et al. 1999; Burrows et al. f997, and 
references therein) . Because of the complex chem- 
istry occurring in their atmospheres, therefore, a 
comprehensive theory of SMOs requires detailed 
knowledge of the opacities of all the abundant 
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molecular species. Obtaining complete opacity 
data has been a major challenge for the field of 
SMO spectral synthesis (Burrows et al. 2001). 

The chemistry of brown dwarf and giant planet 
atmospheres is further complicated by the con- 
densation of gaseous molecules into liquid or solid 
cloud particles, a process which occurs naturally at 
low temperatures. Because clouds significantly af- 
fect spectral features, a satisfactory theory for the 
structure of substellar atmospheres must address 
cloud particle formation and subsequent growth. 
Detailed knowledge of the distribution of particle 
sizes near the photosphere is a basic requirement 
for properly modeling the optical effects of clouds 
(Lunine et al. 1989; Ackerman & Marley 200f ). 

Clouds influence brown dwarf and giant planet 
spectra in several important ways. First, cloud 
formation can deplete the atmosphere of refrac- 
tory elements that become sequestered in con- 
densed form and then rain out from the upper 
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atmosphere. This effect is hypothesized to be im- 
portant for the interpretation of the L to T dwarf 
spectral transition (Burrows & Sharp 1999; Bur- 
rows et al. 2001). Second, clouds near the photo- 
sphere will have the general effect of smoothing 
out prominent spectral features (Jones & Tsuji 
1997). Third, the optical albedo of irradiated ob- 
jects will be increased substantially by the pres- 
ence of clouds (Sudarsky et al. 2000). Fourth, the 
presence of an optically thick cloud layer causes a 
back- warming effect that results in heating of the 
atmosphere. 

Clouds may be manifest in brown dwarf spec- 
tra in a more subtle way. Bailer- Jones & Mundt 
(2001) reported photometric I-band variability of 
up to 7% in L dwarf spectra. The effect is a possi- 
ble signature of variations due to the patchiness of 
clouds. Unfortunately, cloud patchiness in SMOs 
is not well understood because even though much 
is known about the dynamic meteorology of the 
Earth, these gaseous bodies are dynamically very 
different. For example, general circulation models 
of the Earth, although they are of great utility for 
studying Mars, completely fail to reproduce the 
observed dynamics in Jupiter's atmosphere. A de- 
tailed, three-dimensional cloud model paralleling 
the state of the art in Earth climate models would 
require detailed meteorological measurements for 
proper parameter calibration that are simply not 
available for any other planet. 

The present cloud model at its core employs 
the basic microphysical timescale arguments of 
Rossow (1978) to determine the most probable 
particle size in a cloud at each vertical pressure 
level. The model is one-dimensional only. We 
make no attempt to treat the intricacies of cloud 
patchiness or the effects of winds and horizontal 
advection. The present cloud model is not in- 
tended to calculate the detailed meteorology of 
planetary atmospheres but to offer a simple pre- 
scription for estimating the opacity of clouds in 
SMO atmospheres. A detailed but computation- 
ally cumbersome SMO climate model, if it could 
be developed, would be difficult to incorporate 
into spectral synthesis models. The philosophy be- 
hind a simple, one-dimensional approach to cloud 
modeling is to develop a prescription based on 
well-established physical principles that can be 
used to guide and inform spectral synthesis cal- 
culations. The present model aims to characterize 



SMO cloud particle sizes and densities in a global 
and time averaged sense. We therefore assume the 
particle number density of the cloud to be uniform 
across the spherical surface of an object at a given 
pressure and temperature, and that there has been 
sufficient time in these systems to establish chem- 
ical equilibrium throughout the atmosphere. 

Our timescale arguments compute an average 
effectiveness for each of the various competing 
physical processes. Therefore, the output parti- 
cle sizes clearly depend directly on our assump- 
tions of the values of the four unknown parame- 
ters in the microphysical timescales (see Section 
2.3). The calculation represents a first-order esti- 
mate of SMO cloud particle sizes, and we therefore 
present our results for SMO particle sizes with the 
realization that the procedure employed could po- 
tentially overestimate or underestimate the true 
particle sizes by some unknown factor, which we 
expect will be of order unity. Nevertheless, since 
we have been consistent in the assumed values of 
the unknown physical parameters throughout the 
calculation, we expect that the errors in the com- 
puted particle sizes will be uniform among the 
SMOs studied. We are therefore confident, de- 
spite the inaccuracies of the approach, that the 
sizes computed are correct to within an order of 
magnitude and that the trends we demonstrate, in 
which typical cloud particle sizes vary systemati- 
cally with brown dwarf effective temperature and 
gravity, do represent physically meaningful results. 

In this paper, we reiterate the general conclu- 
sion arrived at in Luninc et al. (1989) and Acker- 
man & Marley (2001): it is not satisfactory when 
modeling brown dwarf spectra to assume a pri- 
ori a single particle size distribution because the 
sizes of cloud particles vary strongly with effec- 
tive temperature, surface gravity, and height in 
the atmosphere. We extend the previous efforts 
to incorporate clouds into spectral models by cal- 
culating particle sizes based on a one-dimensional 
model of cloud particle growth for several abun- 
dant species over a range of realistic atmospheric 
conditions. We compute particle radii spanning 
a broad range from about 5 /zm to over 300 /im. 
Therefore, in the context of this discussion, we 
hereafter refer to particles less than 10 /xm in ra- 
dius as small; we refer to particles in the range 
from 10 — 100 /jm in radius as medium-sized; and 
we consider large particles to be those having radii 
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greater than 100 /im. 

In Section 2, we present an improved model of 
cloud formation and droplet growth to determine 
the composition, abundance, and distribution of 
cloud particles in brown dwarf atmospheres. In 
Section 3, we obtain modal cloud particle sizes 
for several representative cloud-forming species 
of high abundance for atmospheric models span- 
ning a broad range of effective temperatures and 
surface gravities. In Section 4, we compare our 
particle sizes with those of two other recent pa- 
pers that have also addressed clouds (Ackerman 
& Marley 2001; Helling ct al. 2001). In Section 5, 
we apply our cloud model to brown dwarf spec- 
tra by computing the opacity of clouds result- 
ing from the wavelength-dependent effects of Mie 
scattering and absorption of radiation. We also 
demonstrate the importance of clouds as an opac- 
ity source by comparing the spectrum of a cloudy 
brown dwarf atmosphere with that of a cloud-free 
atmosphere. 

2. The Cloud Model 
2.1. Condensation Level 

The present cloud model follows much of the 
formalism of Lewis (1969), Rossow (1978), Steven- 
son & Luninc (1988), and Luninc et al. (1989). To 
aid the reader, we list in Table 1 all the symbols 
used in our cloud model as described in this sec- 
tion of the text. 

As in Lewis (1969), we compare the partial 
pressure of a given condensable species to its sat- 
uration vapor pressure. Only at levels where the 
partial pressure exceeds the saturation vapor pres- 
sure can condensation begin. However, the re- 
quired level of supersaturation for efficient con- 
densation depends on the physics of the nucleation 
process itself. We discuss the two relevant nucle- 
ation processes — homogeneous and heterogeneous 
nucleation — in more detail in Section 2.2. In a 
clean atmosphere, in which the availability of con- 
densation nuclei is low (the case of homogeneous 
nucleation), condensation will not begin until the 
vapor becomes highly supersaturated. Thus, the 
partial pressure will greatly exceed — often by a 
factor of two or more — the saturation vapor pres- 
sure (Rogers & Yau 1989). This situation is hy- 
pothesized to occur in methane clouds in Titan's 
atmosphere (Tokano et al. 2001). The condensa- 



tion level appears where 

Pcond > Psat (1 + »$max)- (1) 

In Equation 1, P con d is the partial pressure of 
the condensing vapor, and P sa t is its saturation 
vapor pressure. The <S max parameter is the as- 
sumed maximum supersaturation. For example, 
if <S max = 0.01, then the partial pressure of the 
condensing vapor must be 1% larger than the satu- 
ration vapor pressure before condensation can be- 
gin. To aid future discussions, we also define the 
saturation ratio as S = P CO nd/Psat- 

In addition to the original Lewis (1969) cloud 
model, we also treat the case of heterogeneous (or 
chemical) condensation, in which chemical con- 
stituents different in composition from the cloud 
react chemically to produce the condensate. For 
example, we allow Mg, Si, and O to produce 
forsterite, Mg 2 SiC>4. In such cases of heteroge- 
neous condensation, the concept of a saturation 
vapor pressure is replaced by the heterogeneous 
condensation pressure. We apply the chemical 
equilibrium model of Burrows & Sharp (1999) to a 
solar-composition gas (Anders & Grevesse 1989) 
to determine which chemical species are thermo- 
dynamically favored. The calculation determines 
the equilibrium compositions by minimizing the 
Gibbs free energy of the system. 

The maximum supersaturation, iS max in Equa- 
tion 1 , appears in the calculation as an essentially 
free parameter, though we have some guidance 
as to its likely value from measurements of the 
supersaturation of water clouds on Earth. We 
assume in these calculations that the abundance 
of condensation nuclei is sufficient for heteroge- 
neous nucleation to begin at only 1% maximum 
supersaturation, though we will explore the ef- 
fect of this parameter on the computed particle 
sizes subsequently (see Section 3.4). We have not 
attempted herein to calculate the precise com- 
position or abundance of these condensation nu- 
clei. But even for iron, which is the most re- 
fractory of the homogeneously condensing species, 
we argue based on chemical equilibrium considera- 
tions that there exist abundant refractory conden- 
sates that can potentially serve as nucleation sites 
for heterogeneous nucleation, including corundum 
(AI2O3), hibonite (CaAli 2 Oig), the calcium ti- 
tanates (e.g., CaTiC>3 and Ca4Ti 3 Oio), grossite 
(CaAUOy) and other calcium-aluminum oxides, 
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Table 1 

Cloud Model Variables and Parameters 



Symbol Description 



g sur f, T e ff SMO surface gravity and effective temperature. 

Pcond, Psat Partial pressure of the condensing vapor and its saturation vapor pressure. 

Smax, S Assumed maximum supersaturation and the saturation ratio. 

Tcond, T nuc Timescales for heterogeneous and homogeneous nucleation. 

Tcoal; T coag Timescales for coalescence and coagulation. 

Tail, T conv Timescales for gravitational fallout and convective upwelling. 

r ra d Timescale for particles to cool by radiation. 

T, P, p, i] Temperature, pressure, mass density, and dynamic viscosity of the atmosphere. 

R, kb,NA Universal Gas Constant, Boltzmann's Constant, and Avogadro's Number. 

H, /ip Mean molecular weight of atmosphere and condensate molecular weight. 

Kn, Re Knudsen and Reynolds numbers. 

A, c soun d Mean free path of atmospheric gas particles and atmospheric sound speed. 

v Relative velocity between particle and fluid in Re expression. 

Vfaii(r), g Terminal velocity of particles of radius r in local acceleration of gravity, g. 

Vconv Upward velocity of convecting gas parcels. 

H Atmospheric scale height; i.e., pressure e- folding distance. 

e sur f , L Surface tension of condensed molecules and the latent heat of vaporization. 

r c , Z Critical radius of particles and the Zcldovich factor of classical nucleation theory. 

N Number density of cloud particles. 

£coai ; fcoag Coalescence and coagulation efficiencies. 

q c Moles of condensate per mole of atmosphere. 

q v Moles of condensable vapor per mole of atmosphere. 

q t Total moles of condensable material per mole of atmosphere. 

qbciow Mole fraction of condensable vapor below the cloud base. 

P c i Condensation curves as shown in Figure 1. 

n(r), ro Particle size distribution about modal particle size, ro. 
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and spinel (MgAl 2 4 ). 

As in Rossow (1978), we ignore variations in 
<S ma x within the cloud. Although in reality the 
supersaturation should decrease with altitude, we 
ignore this effect in the present calculations. Fur- 
thermore, it is possible to have a small amount 
of condensate present between S = 1 and S = 
1 +>5 max , even though condensate cannot possibly 
form there (since the initial condensation of vapor 
into embryonic particles in general requires at least 
a modest level of supersaturation) . However, par- 
ticles can settle into the region between S = 1 and 
S = 1 + iS max once they have formed and yet still 
be thermodynamically stable. Since it depends on 
a balance between the processes of gravitational 
settling and convective upwelling, the amount of 
cloud material present in the barely supersatu- 
rated region is difficult to calculate with certainty. 
The clouds we treat in this article have a very low 
<Smax, and therefore the errors introduced by ne- 
glecting these effects are small. In a cloud having 
very few condensation nuclei, which as we have 
stated would require a large <S max , considerations 
such as the variation of the supersaturation with 
altitude and the exact location of the cloud base 
become more crucial for accurate cloud modeling. 
The present cloud model is not configured to ac- 
count for these effects realistically, though they 
are relatively unimportant for the low supersatu- 
rations generally assumed throughout this paper. 

Figure 1 shows the condensation curve for two 
condensates, forsterite (Mg2Si04) and gchlenite 
(Ca2Al2SiC>7), as well as the "effective condensa- 
tion curve" of iron derived from the iron saturation 
vapor pressure curve, and several model brown 
dwarf atmospheres. As the caption to Figure 1 
explains, the curve labeled "Fe" is not the true 
saturation vapor pressure curve of iron. It is an 
"effective condensation curve" derived by taking 
the saturation vapor pressure of iron and dividing 
by the number mixing ratio of iron vapor just be- 
low the cloud base (parameter qbeiow of Table 2), 
which is the maximum allowed number mixing ra- 
tio of gaseous iron in a solar composition mixture 
when hydrogen is all in molecular form (H2). We 
discuss this further in Section 2.3. 

The intersection points between the four long- 
dashed curves of brown dwarf temperature- 
pressure profiles and the condensation curves of 
forsterite and gehlenite represent levels in the at- 
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Fig. 1. — Illustrates the locations of the cloud 
bases for four different atmospheric models (long- 
dashed curves) and three different chemical species 
(solid and short-dash curves). The atmospheric 
models are at two different surface gravities and 
two different effective temperatures. The short- 
dashed curve defines the atmospheric pressures 
and temperatures for which the vapor pressure 
of iron equals its saturation vapor pressure. This 
curve is obtained using the iron saturation vapor 
pressure expression in Barshay & Lewis (1976). 
The saturation vapor pressure curve has been 
adjusted using the iron vapor mixing ratio (by 
number) below the cloud (parameter qbeiow in 
Table 2). The two solid lines arc the curves 
which illustrate when condensation of forsterite 
(Mg 2 SiC>4) and gehlenite (Ca2Al2SiC>7) become 
thermodynamically favored. Thus, they repre- 
sent the pressure dependent condensation temper- 
atures of those species under the assumption that 
chemical equilibrium holds. The intersection of a 
given atmospheric model and the two condensa- 
tion curves represents the pressures and temper- 
atures at the bases of the chemically condensing 
clouds in our cloud model. However, for iron, the 
cloud base will be at lower temperature and pres- 
sure than the intersection point if the iron vapor 
is significantly supersaturated (see Equation 1). 
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mosphcrc at which the relevant species may ap- 
pear via one or more chemical reactions. 

The more familiar case of homogeneous con- 
densation, or direct condensation from a vapor of 
the same composition, does occur for water and 
iron. Thus, for the condensation of iron and wa- 
ter, we follow Lewis (1969) and employ the well- 
known saturation vapor pressure relations of iron 
and water (Barshay & Lewis 1976; Lunine et al. 
1989). In Figure 1, the iron curve is thus shown as 
a thick, short-dashed line as a reminder that the 
iron cloud condenses directly from iron vapor. As 
we discuss in more detail in Section 5.3, in which 
we explore the coupling between clouds and the 
radiative transfer problem, the atmospheres used 
in this stage are computed on the basis of an in- 
dependent, dust-free stellar atmosphere code. 

2.2. Microphysical Timescales 

Although the chemistry applying to brown 
dwarf atmospheres is radically different from 
Earth's atmospheric chemistry, we nevertheless 
treat the various competing effects in the cloud 
formation process as microphysical timescales 
(Rossow 1978). In our model, five microphysi- 
cal processes compete with gravitational fallout 
to increase the modal particle size: 1) convec- 
tive uplifting of condensable vapor, 2) heteroge- 
neous nucleation, 3) homogeneous nucleation, 4) 
coagulation, and 5) coalescence. Each process is 
characterized by a single-valued timescale which 
expresses its relative importance. 

The expressions for four of the microphysical 
timescales are computed in Rossow (1978): r con d 
(heterogeneous condensation), r coa i (coalescence), 
Tcoag (coagulation), and Tf a ii (gravitational fall- 
out). These expressions depend on the values of 
both the Knudscn and Reynolds numbers, which 
characterize how particles of different sizes inter- 
act physically with the surrounding fluid and with 
each other. It is therefore necessary to calculate 
these values during the particle growth phase of 
the calculation (see Section 2.3). 

We also require an equation of state relating 
the gas density to its temperature and pressure. 
In this model, we employ the ideal gas equa- 
tion of state to relate the gas mass density p to 
the atmospheric temperature and pressure (T, P): 
p = /iP/RT, where R is the universal gas constant, 



8.314 x 10 7 ergsmol _1 K _1 . The equation of state 
depends on the mean molecular weight of the at- 
mospheric gas mixture. For these atmospheres, 
we use p = 2.36 gmol -1 , the mean molecular 
weight of a solar composition gas in which hydro- 
gen is present in the molecular state, H 2 . We have 
used the chemical equilibrium model of Burrows 
& Sharp (1999) to verify that — as for Jupiter — 
molecular hydrogen is the dominant phase of hy- 
drogen in the upper atmospheres of SMOs. The 
ideal gas equation of state reproduces to within 1% 
the value of p obtained using the more sophisti- 
cated thermodynamic equation of state described 
in Burrows et al. (1989). 

The dimensionless Knudsen number differenti- 
ates between the two regimes — classical and gas 
kinetic — of physical interaction between cloud par- 
ticles and the surrounding medium: 

Kn = A/r, (2) 

where r is the cloud particle's radius and A is the 
mean free path of atmospheric gas particles. The 
transition between the two regimes occurs where 
the Knudsen number is near unity (i.e., where 
A ~ r). The high Knudscn number regime, in 
which the mean free path of gas molecules is much 
larger than the size of a typical cloud particle, 
will normally not occur in SMO clouds. However, 
there are two situations in which large Knudscn 
numbers could arise: 1) when particles are smaller 
than ~1 pm in radius, and 2) high up in the at- 
mosphere of a planet or brown dwarf at pressures 
much less a bar. We include the Knudsen number 
in the cloud model to account for these possibili- 
ties. 

The dimensionless Reynolds number character- 
izes the behavior of a spherical particle moving 
through a fluid. The equation for the Reynolds 
number is 



where p is the mass density of gas in the atmo- 
sphere, r is the radius of the particle, v is the 
velocity of the particle relative to the fluid, and 
r\ is the dynamic viscosity of the fluid. The flow 
will be turbulent for Re 3> 10 but laminar for low 
Reynolds numbers (Landau & Lifshitz 1959). For 
the purposes of these calculations, we assume lam- 
inar flow unless Re > 70 (see Section 2.3). This 



6 



choice for the laminar to turbulent cutoff is con- 
venient because above Re = 70, the drag coeffi- 
cient describing the variation in terminal velocity 
of a falling particle becomes roughly a constant at 
C D « 0.2 (Rossow 1978). 

To calculate the Reynold's number and the 
timescale for gravitational fallout, we introduce 
an approximate value for r/, the dynamic viscos- 
ity. The viscosity of a hydrogen and helium gas 
mixture at low densities is roughly constant over 
a broad range of temperatures. Because the vis- 
cosity increases with the square root of temper- 
ature at low pressures, it varies by less than a 
factor of two over the full range of SMO effec- 
tive temperatures discussed herein (see Section 
2.3). For simplicity, our model therefore assumes 
a constant value for the dynamic viscosity of rj = 
2 • 10 -4 poise, which is consistent with the mea- 
sured values of 77 for hydrogen and helium at low 
densities and T - 1000 K (Lide 1991, CRC Hand- 
book). 

The Tf a ii timescale represents gravitational fall- 
out of particles, which is the fundamental process 
limiting particle growth. As the particles become 
large, they fall out. We take the gravitational fall- 
out timescale to be the time required for particles 
to fall through a pressure scale height of the at- 
mosphere at terminal velocity. 

The terminal velocity of falling particles de- 
pends on their size, shape, density, rigidity, and 
the nature of the fluid flow around them (e.g., 
Stokes flow, free molecular flow, or turbulent flow). 
We assume rigid spherical particles throughout 
this paper. While it is true that liquid particles 
of a given size will fall more slowly through the 
atmosphere than solid particles because of the dis- 
ruption of their shape, the effect will not signifi- 
cantly alter the fall speeds for particles smaller 
than about a millimeter (Rogers & Yau 1989). 
The particles we obtain using this model — see Sec- 
tion 3 — never grow large enough for the difference 
in terminal velocity between solid and liquid par- 
ticles to be noticeable. Thus, taking the particles 
to be rigid spheres is in this case a reasonable sim- 
plification. 

For the common case of low Knudsen and low 
Reynolds number flow, the terminal velocity will 
be given by the Stokes solution for a falling rigid 
sphere (Rogers & Yau 1989). The Cunningham 
factor, 1 + f Kn, accounts for the variation in the 



terminal velocity with Knudsen number. There is 
also a variation in terminal velocity with Reynold's 
number. These combined effects yield three differ- 
ent forms for the terminal velocity Vf a n as a func- 
tion of the radius r of the particle: 

v fa „(r) = r 2 . (^f^ 1 ) Rc<70, Kn < 1, 

j V __ (4) 

vfaii(r) = r 1/2 y 40P 3 C °' ldg Rc > 70, Kn < 1, and 

, P _ (5) 

VfaU « =r 7^^b^J Kn>L 

(6) 

In Equations 4, 5, and 6, Na is Avogadro's 
number (6.02 2 x 10 23 mol _1 ), kb is Boltzmann's 
constant (1.38 x 10 _16 crgK _1 ), and g is the local 
gravitational acceleration. The quantity p CO nd rep- 
resents the mass density of the particle itself, for 
which we use the bulk density of a solid or liquid 
phase having the composition of the condensate. 
The parameters used for each condensate we have 
treated — iron, forsterite, gehlenite, and water — 
are listed in Table 2. The mass density of the gas p 
depends on temperature and pressure through the 
equation of state. The terminal velocity of parti- 
cles depends also on the local acceleration of grav- 
ity, g, which varies considerably among the SMO 
population (see Section 3). In the Stokes solution, 
Equation 4, the terminal velocity is proportional 
to the square of the radius of the falling sphere. 
In the turbulent regime, Equation 5, it is propor- 
tional to the square root of the radius. In the gas 
kinetic (or large Knudsen number) regime, Equa- 
tion 6, the terminal velocity increases linearly with 
particle radius. The gravitational fallout timescale 
is inversely proportional to the terminal velocity: 
Tfaii = H/vf a n, where H = ^ is the atmospheric 
scale height. We describe in Section 2.3 our pro- 
cedure for differentiating between the Reynolds 
number regimes, which is necessary to calculate 
Tfaii and the other microphysical timescales. 

The T conv timescale represents the convective 
upwelling of particles. Although convective up- 
welling is not itself a particle growth process, con- 
vective updrafts hold small cloud particles aloft in 
the atmosphere so that they accumulate more ma- 
terial via coagulation, coalescence, etc. We define 
the timescale characterizing convective updrafts as 
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Tconv = H/v conv . The convective timescale is thus 
the time required for gas parcels to rise through 
a pressure scale height. The upward convective 
velocity, v conv , is estimated via a mixing length 
prescription in the atmosphere code. 

The nucleation timescale, r nuc , characterizes 
particle growth via homogeneous nucleation, in 
which highly supersaturated vapor condenses 
spontaneously due to rapid collisions between the 
vapor molecules. For the nucleation timescale, 
we employ the classical homogeneous nucleation 
theory, in which droplets form in free space as a 
result of chance collisions between the molecules 
of the supersaturated vapor. This process is rarely 
observed to occur in nature because its timescale 
will be longer in almost all cases than r con d , which 
characterizes the condensation growth of a pop- 
ulation of particles via heterogeneous nucleation. 
Condensation growth means the conversion of va- 
por to particles, whether they be in the liquid 
or solid phase. For example, the phase diagram 
of iron indicates that iron will condense into liq- 
uid particles above 1800 K. However, we ignore 
the distinction between liquid and solid particles 
because the particles are small enough that the 
flattening effect during free-fall is negligible. In 
heterogeneous nucleation, particles nucleate onto 
seed particles, called condensation nuclei, which 
are presumed to populate the region of conden- 
sation. As discussed, when heterogeneous nucle- 
ation is efficient, condensation may begin even 
when the level of supersaturation is quite low 
(< 1%). Neither nucleation process represented 
by T nuc and r con d apply to the species that ap- 
pear due to chemical reactions (e.g., forsterite and 
gehlenite) rather than from a supersaturated va- 
por. The nucleation timescale expressions apply 
only to the species (e.g., iron and water) for which 
the condensate forms directly from a vapor of the 
same composition. 

For the condensation timescale, T conc j , we con- 
sider nucleation onto seed aerosols, also known 
as condensation nuclei. The assumption of a 
low maximum supersaturation of the condensable 
vapor — <S max = 1% throughout (see Section 2.1) — 
is equivalent to the assumption that condensation 
nuclei are abundant in the upper atmosphere; i.e., 
these are "dirty" atmospheres. We have not at- 
tempted in the present model to calculate the com- 
positions or exact abundances of these condensa- 



tion nuclei. Our assumed value for <S max is based 
rather on the recognition that the complex chem- 
istry occurring in SMO atmospheres will likely 
generate a plethora of molecules suitable as seeds 
for the onset of heterogeneous nucleation. The ex- 
pressions for the condensation timescale, which de- 
pend on the Knudscn number regime, are adapted 
from Rossow (1978): 



Tcond = pr 



PcondRT 
47?M P sat'5 max 



Kn < 1, (7) 



Tcond 



/ PcondRT 
\ 3Psat^max 



2/xN A k b T 



1/2 



Kn > 1. 

(8) 



In Equations 7 and 8, P sat is the saturation va- 
por pressure defined in Equation 1, and p con d is 
the mass density of the condensate. The other 
variables are the same as in Equations 4, 5, and 6. 
The mass densities of the four condensates treated 
in this paper are given in Table 1. 

The condensation timescale given in Equa- 
tions 7 and 8 assumes that heterogeneous nu- 
cleation dominates over homogeneous nucleation; 
i.e. c> max <C 1. We relax this assumption by incor- 
porating the process of homogeneous nucleation 
explicitly (r nuc ). It should be noted, however, that 
the inclusion of T nuc will only be important if the 
maximum supersaturation parameter is assumed 
to be large (<S max S> 1). We present homogeneous 
nucleation as a feature of the cloud model so that 
particle growth can be treated in the extreme case 
of very high supersaturations. As we have indi- 
cated, however, this situation is not likely in brown 
dwarf atmospheres. In most cases, homogeneous 
nucleation will be orders of magnitude slower than 
condensation growth by heterogeneous nucleation 
(i.e., r nuc <C T cond ). 

For the timescale of homogeneous nucleation, 
we have adapted expressions from Stevenson & 
Lunine (1988) for large Knudsen numbers: 



^"nuc 



PcondrL 2 ^ 



[exp (2e sur f^/RTp cond r) - 1] 



(Kn> 1). 
(9) 



Here, r is the particle radius, L is the latent heat 
of vaporization, e sur f is the surface tension of the 
condensed liquid molecules, c soun( j = ^/RT//x is 
the sound speed in the atmosphere, and P and T 
are the ambient temperature and pressure. We list 
the values of these parameters for iron and water 
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in Table 2. This expression is derived from Equa- 
tions 3-5 of Stevenson & Lunine (1988) by tak- 
ing r nuc = r and then substituting in from 
the ideal gas equation of state and the Clausius- 
Clapeyron equation to eliminate the particle num- 
ber density, saturation vapor pressure, and the 
vapor pressure temperature gradient. We take 
the efficiency of heat exchange in Equation 3 of 
Stevenson & Lunine (1988) to be 100% and the 
average relative velocity of cloud particles collid- 
ing with the local hydrogen to be c soun d- 

The T nuc expression for large Knudscn num- 
bers represents the timescale for the growth of 
an embryonic cloud particle that must release 
the latent heat of condensation through collisions 
with the surrounding hydrogen (Stevenson & Lu- 
nine 1988). It is possible, however, for parti- 
cles in high-temperature environments to grow by 
releasing the latent heat of condensation simply 
through radiatively cooling. We therefore include 
in the model an estimate of the radiative cooling 
timescale, which is compared with Equation 9 in 



-l 

'rad 



r-1 



path of particles is greater than the particle size, 
we use T ra( j in place of r nuc . The equation we use to 
approximate r ra d is adapted from Woitke (1999): 



(2 x 10" 



sec 



Pcond 



1 



1 gmol 1 



; cm" 
T 



1000 K 



(10) 



In Equation 10, /i p is the molecular weight of the 
condensate (as distinct from fj, w 2.36 gmol -1 , 
the mean molecular weight of the atmospheric gas 
mixture). The values of /i p for the four conden- 
sates treated herein are given in Table 2. 

It should be noted that in a very high temper- 
ature, low density environment, in which the re- 
lease of latent heat is extremely fast, the limiting 
timescale for homogeneous nucleation will be nei- 
ther T nuc nor r ra( j but the timescale for the diffu- 
sion of vapor molecules to the surface of the grain. 
We do not calculate this timescale in the present 
model. For the purposes of SMO atmospheres, 
however, this regime will never be realized in 
practice because the condensation temperatures of 
even the most refractory compounds rarely exceed 
2000 K. Therefore, although the limit of r ra d — ► 
in this formulation appears to be potentially prob- 



lematic, the temperature constraints ensure that 
r ra d cannot be arbitrarily small. Hence, for use 
in the calculation of SMO clouds, it is reasonble 
to assume that the shorter of the two timescales, 
r nuc and T ra d, determines the particle's homoge- 
neous nucleation rate. 

The constant preceding Equation 10 is strictly 
correct only for highly opaque species having high 
extinction efficiencies. This is appropriate for 
the condensates we focus on here, although the 
timescale for radiative cooling will be longer by 
a factor of ten or a hundred than the timescale 
shown in Equation 10 if the particles are mostly 
transparent. In cases in which the T ra d 

timescale 

dominates the particle sizes, a more accurate ap- 
proximation than we have made in Equation 10 
should be calculated by integrating the extinction 
efficiency over wavelength; e.g., see Equation 17 of 
Woitke (1999). 

The more common case of small Knudsen num- 
bers, which corresponds to high atmospheric gas 
density, is a commonly used meteorological expres- 
sion; see Rogers & Yau (1989). The low Knudscn 
number expression for T nuc is written in terms of 
the critical radius at which the equilibrium vapor 
pressure over the surface of a particle equals the 
ambient vapor pressure. The equilibrium vapor 
pressure over the surface of a liquid particle, be- 
cause of its finite curvature, is in general higher 
than the equilibrium vapor pressure over a flat 
surface of the liquid, which is the saturation va- 
por pressure normally measured in the laboratory. 
Therefore, the critical radius of the liquid parti- 
cle will only be attained when the saturation ratio 
exceeds unity; i.e., when the actual ambient pres- 
sure exceeds the saturation vapor pressure. The 
critical radius therefore depends on the supersat- 
uration 5 max (Rogers & Yau 1989): 



2/^p€surf 



nd RTln(l + 5 n 



(11) 



Unlike the gas kinetic regime, in which the re- 
moval of latent heat limits the rate at which par- 
ticles can grow by nucleation, in the classical gas 
regime, the nucleation rate of particles is deter- 
mined by the rate at which supercritical droplets 
are formed; i.e., droplets larger than the critical 
radius above which particles grow spontaneously. 
The expression for T nuc in terms of the critical 
droplet radius r c , adapted from Rogers & Yau 
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(1989), reads 



of the condensed particle, r, is given by 



— T L exp 



87rN A r4 



47rrge surf 
3k b T 



N 



(Kn < 1) 
(12) 

In Equation 12, Z is the dimensionless Zcldovich 
or non equilibrium factor, which depends on tem- 
perature and the physical properties of the con- 
densate. The equation for Z, given by Jacobson 
(1999), is 



27rr2 /9cond N A V k bT 



e surf 



(13) 



The Zeldovich factor accounts for the differences 
between equilibrium and non-equilibrium cluster 
concentrations in the classical homogeneous nucle- 
ation theory used to derive Equations 11 and 12. 

Coagulation and coalescence represent parti- 
cle growth due to collisions. The coagulation 
timescale, r coag , refers to the formation of larger 
particles by the collision of smaller particles. It 
therefore depends primarily on the thermal tem- 
perature, viscosity, and the number density of 
cloud particles. The coalescence timescale, T coa i, 
characterizes the growth of particles by coales- 
cence, in which large particles having high down- 
ward velocities in the fluid overtake and merge 
with small, slowly falling particles. Thus, coales- 
cence refers to the collisional process caused by 
the different fall speeds of different sized particles, 
whereas coagulation is collisional growth resulting 
from Brownian motion. Both processes proceed at 
a rate proportional to the particle number density, 
N, since the total amount of material available for 
condensation is conserved. The number density of 
cloud particles is given in terms of the gas density 
as 

p\ ( ^p 



N 



q c - 



(14) 



V47r/9 con dr 3 

In Equation 14, r is the cloud particle radius, p is 
the mass density of the surrounding gas (as given 
by the equation of state) , p con d is the mass density 
of the condensate, and q c is the condensate mixing 
ratio by number, which depends on the specific 
rainout prescription applied. The rainout scheme 
we employ to compute the variation of q c with 
height above the cloud base is explained in Section 
2.3. 

Given the particle number density, N, the 
timescale for coagulation in terms of the radius 



' coag 



dN/dt 



3idL;w' Kll<1 (15 » 



' coag 



4e r 



Pcond 1 Kn > L 



ccoag v 3rk b T N' " (16) 

Here, e coag is the coagulation efficiency, which is 
a free parameter of the model. We discuss this 
further in Section 2.4. 

Our final microphysical timescale, T coa i, charac- 
terizes the process of coalescence, in which larger 
particles overtake and coalesce with smaller par- 
ticles. For coalescence, the Reynolds number 
again plays an important role because — unlike 
coagulation — coalescence pertains specifically to 
the merging of particles of greatly different radii 
whose interaction depends on the nature of the 
surrounding flow fields. The coalescence expres- 
sions are adapted from Rossow (1978), although 
we have introduced an extra free parameter into 
the equations, the efficiency for coalescence, e coa i: 



1 



Tc ° al r 4 e coal V7rp cond gy N' 



9t7 



1 



Re < 70, Kn < 1 



7"coal 



7"coal 



3p 



£coal V 107T 2 p cond r 5 g N 

1 / 27/j \ [ 2RT I 



(17) 

Re > 70, Kn < 1 
(18) 



I" ^coal 



Kn > 1. 



(19) 

As we discuss in Section 2.4, coalescence is im- 
portant only for liquid particles. This is the only 
case in our model in which the distinction between 
liquid and solid particle formation is important. 
Thus, for iron cloud formation above 1800 K, coa- 
lescence will potentially be important, and it can 
be important for water clouds, but we ignore it for 
the silicate clouds. 

2.3. Our Cloud Code 

The code proceeds in two stages: (1) determina- 
tion of the cloud base globally in the atmosphere, 
and then (2) particle growth at each atmospheric 
temperature-pressure level. The algorithm is run 
for each of the atmospheric temperature-pressure 
profiles, which differ in effective temperature and 
gravity, and for each condensable species sepa- 
rately. In our model, we ignore core-mantle grains 
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composed of molecules of two or more of the major 
species. These types of condensates could poten- 
tially be important in regions in which several ma- 
jor constituents in the equilibrium vapor mixture 
form at similar temperatures and pressures. In the 
present cloud model, we do not treat the growth 
of grains composed of multiple chemical phases. 

We list in Table 2 the input parameters for the 
four condensable species treated herein. 

In stage (1), two cases must be treated sepa- 
rately as explained above. In the case of homo- 
geneous condensation, we calculate the intersec- 
tion between the partial pressure of the condens- 
able vapor and its saturation vapor pressure curve 
and identify the condensation level using Equation 
1. In the case of heterogeneous condensation, in 
which the condensing chemical is not present in 
vapor form below the cloud base, we determine 
the cloud base by using the chemical equilibrium 
model of Burrows & Sharp (1999) (see Figure 1). 

In stage (2), we calculate the modal size of 
particles, which are assumed to be spherical, at 
each atmospheric temperature-pressure level. The 
model outputs the modal particle radius. This 
radius is determined by growing embryonic par- 
ticles of radius ~10A slowly until the various 
competing particle growth timescales are balanced 
by the timescale of gravitational fallout. For 
the most probable particle radius, the sedimen- 
tation timescale equals the shortest of the growth 
timescales: 

Tfall = Min {T nuc , T con d, T coa g, T coa l, T conv }. (20) 

Note that for homogeneously condensing species 
in the gas kinetic regime (Kn > 1), it may be that 
radiative cooling, rather than the release of latent 
heat by collisions with atmospheric hydrogen, is 
the dominant process limiting the rate of parti- 
cle growth. Therefore, if the Knudscn number is 
greater than one and r r ~ d > t~* c , we replace r nU c 
in Equation 20 with the timescale for particles to 
release the latent heat of condensation by radia- 
tion, T rac i. Additionally, as we have stated in Sec- 
tion 2.2, nucleation does not apply to the case of 
chemical clouds (Rossow 1978), and coalescence is 
inefficient between solid particles. The condensed 
particles can only accumulate more material by 
coagulating under the influence of convective up- 
drafts. Hence, for forstcrite and gehlenite, the fall 
timescale must balance the minimum of the two 



timescales relevant for heterogeneously condens- 
ing clouds: T coag and r conv . 

The timescales in Equation 20 depend on the 
Knudsen and Reynolds numbers (see Section 2.2). 
We therefore need a prescription to determine 
which physical regime applies as the particles are 
grown. Because the Reynolds number depends on 
velocity, and the velocity depends on the Reynolds 
number, calculating the two quantities indepen- 
dently for a given particle size and atmospheric 
parameters is a circular problem. There arc no 
difficulties if the Knudsen number is larger than 
one, since the gas kinetic terminal velocity in that 
case does not depend on the Reynolds number; the 
velocity to use in this case is given by Equation 6. 
To differentiate between Reynolds number regimes 
in the small Knudsen number case (i.e., Equations 
4 and 5), we use the following procedure based on 
the fact that each of the two expressions for Vf a n 
can be calculated independently of the Reynolds 
number. 

At each point in the particle growth phase, we 
use both Equations 4 and 5, along with the known 
atmospheric parameters and the current particle 
size. In solving for the terminal velocity, we sub- 
stitute the SMO surface gravity, g SU rf, in for g, 
the local gravitation acceleration. We then calcu- 
late two temporary values for the Reynolds num- 
bers, one for each of the velocities calculated. If 
these two independently computed Reynolds num- 
bers are both smaller than our Reynolds number 
cutoff point of 70 between laminar and turbulent 
flow (Rossow 1978), we take the flow around the 
particles to be laminar and adopt the Stokes ter- 
minal velocity, Equation 4, as the actual terminal 
velocity. Likewise, if they are both larger than the 
cutoff, we take the flow around the particles to be 
turbulent and adopt the turbulent terminal veloc- 
ity, Equation 5, as the actual terminal velocity. 

The intermediate case, in which one temporary 
Reynolds number is larger than the cutoff but the 
other is smaller than the cutoff, is more difficult 
to resolve. We settled on the choice of using the 
lesser of the two terminal velocities as the actual 
terminal velocity. Although this could have been 
handled in a variety of ways, the inherent ambigu- 
ities could not have been eliminated without de- 
tailed knowledge of the drag coefficient, which un- 
fortunately depends in general on the Reynolds 
number. Ours was the most conservative choice 
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Table 2 

Condensate Physical Properties 



Model Parameter 


Iron 


Forsterite 


Gehlcnitc 


Water 


Molecular weight, /x p [gmol -1 ] 


55.8 


140.7 


274.2 


18.0 


Mass Density, p con d [gem -3 ] 


7.9 a 


3.2 


3.0 


1.0 a 


Mole Fraction, qbciow b 


5.4 x 1(T 5 


3.2 x 10~ 5 


1.8 x 10~ 6 


1.4 x 10^3 


Surface Tension, e sur f [ergscm~ 2 ] 


200 






75 


Latent Heat of Vaporization, L [ergsg -1 ] 


6.34 x 10 10 






4.87 x 10 12 



a We are ignoring here the minor differences between the densities of the liquid and solid phases 
of iron and water (Lide 1991, CRC Handbook). 

b The qbciow parameter is the mole fraction of condensable vapor just below the cloud base; the 
calculation of qbciow assumes a solar abundance distribution of the elements (Anders & Grevesse 
1989). 



Note. — Values for the surface tension and latent heat of vaporization are taken from Lide (1991, 
CRC Handbook). They are used to derive T nuc using the classical homogeneous nucleation theory 
and are therefore not applicable to the heterogeneously condensing species (see Section 2.2). Note 
also that the mole fractions listed apply only to the cloud base (i.e., the condensation level, as 
shown in Figure 1). Above the cloud base, the mole fraction of total condensable material in the 
cloud drops off much more steeply than the gas pressure (as prescribed in Equations 21 and 22). 
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for a first-order estimate of the terminal velocity. 
With this procedure in place to calculate the ter- 
minal velocity of particles, the effective Reynolds 
number — which is used to calculate all the other 
timescales — is computed directly from its defini- 
tion in Equation 3, with Vf a ii substituted in for 
the relative velocity, v. 

Equation 20 shows the role of convection in this 
model. Convection directly opposes gravitational 
sedimentation. In the presence of sufficiently vig- 
orous convection, therefore, the particles may con- 
tinue to accumulate material until their terminal 
velocities become equal to the convective updraft 
velocity. 

The simple cloud model delineated above is a 
first-order estimate of SMO particle sizes. The 
true particle sizes could very well deviate from the 
values computed by a factor of order unity, though 
we believe that this factor will be the same for all 
objects studied because the systematic errors de- 
pend on our assumptions of the unknown parame- 
ters of the model, assumptions which we have ap- 
plied consistently throughout the computations. 

Our treatment of rainout, in which condensates 
settle gravitationally and thus become depleted 
from the upper atmosphere, differs between homo- 
geneously and heterogeneously condensing species. 
Depletion of the condensate with height must dif- 
fer between these two different types of clouds be- 
cause homogeneous condensation involves a pres- 
sure equilibrium between the vapor and the con- 
densed phase (See Section 2.2), whereas hetero- 
geneous condensation produces the condensate di- 
rectly through exothermic chemical reactions. 

For the purposes of this discussion, we adopt 
the notation of Ackerman & Marley (2001) for 
the relevant mixing ratios: q v = moles of vapor 
per mole of atmosphere, q c = moles of condensate 
per mole of atmosphere. The total mixing ratio is 
qt = q v + q c - We also define qboiow as the mixing 
ratio by number (or mole fraction) of condensable 
vapor just below the cloud base, which equals q t 
at that level, since q c = by definition below the 
cloud base (i.e., at levels in the atmosphere where 
Equation 1 is not satisfied). The values of qboiow 
for each condensable treated herein are given in 
Table 2. 

For both homogeneously and heterogeneously 
condensing species, we assume aggressive rainout 



such that the total material available to condense 
into particles drops off with pressure much more 
steeply than the gas pressure. That is, at the con- 
densation level, the mixing ratio by number (or 
equivalently, the mole fraction) of the condensing 
species is taken to be equal to the maximum par- 
tial pressure attained by that species in a gas of 
solar composition (i.e., parameter qbdow in Table 
2). For heterogeneously condensing species, qbdow 
is constrained by the maximum partial pressure of 
the least abundant of the constituent molecules 
combining to form the condensate. For example, 
in the case of forsterite, the magnesium abundance 
in the gas controls the value of qbdow- We have as- 
sumed a solar abundance of magnesium in the at- 
mosphere below the forsterite condensation level. 
High above the condensation level, however, the 
mixing ratio of cloud-forming material is signif- 
icantly reduced from the equilibrium ratio, and 
the cloud becomes more and more tenuous with 
decreasing pressure. 

For homogeneously condensing species, our 
rainout prescription is based on the assumption 
that all of the supersaturated vapor above the 
base of the cloud will condense. As we have dis- 
cussed in Section 2.1, the onset of particle for- 
mation requires a supersaturated environment, 
which thus elevates the cloud base above the level 
where precise saturation of condensable vapor is 
attained. Once present, however, particles re- 
main thermodynamically stable unless the satu- 
ration ratio drops to below one, which will not 
happen without a temperature inversion in the 
atmosphere's thermal profile. For already existent 
particles or cloud particle embryos, the vapor in 
equilibrium with the particles remains exactly at 
saturation. All excess material goes into the con- 
densate. Hence, throughout the cloud, we take 
qv = Psat/P, where P sat is the saturation vapor 
pressure and P is the gas pressure. Notice that 
this ratio decreases with altitude because P sat de- 
creases with pressure much more steeply than the 
gas pressure, P. It therefore follows that the par- 
tial pressure of the condensable vapor required 
to maintain a saturated environment will be ex- 
tremely small high above the condensation level. 
Though the mole fraction of condensing vapor, 
q v , is constrained by the requirement that the 
condensed phase remain thermodynamically sta- 
ble during particle growth, we still need the total 
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number mixing ratio (i.e., mole fraction) of con- 
densable material, qt, in order to determine q c , 
the condensate number mixing ratio. 

The character of the rainout itself derives from 
the assumption that the total material available 
for condensation above the cloud base follows the 
saturation vapor pressure and not the gas pres- 
sure. We use q t = (l+<S ma x) Psat/P, which implies 
that the total condensable material is constrained 
to drop off with altitude proportionally to the sat- 
uration vapor pressure. This is an assumption of 
the cloud model. It is not based on a rigorous 
derivation of the total mixing ratio with altitude. 
The justification for the rainout scheme is that 
once condensation begins in the atmosphere, with 
subsequent growth and gravitational settling of 
particles, the condensable material can no longer 
remain well-mixed in the atmosphere, as we as- 
sume it had been below the cloud base. There- 
fore, the total number of moles of condensable ma- 
terial per mole of gas will decrease rapidly from 
its initial value with increasing height above the 
condensation level. By contrast with our rainout 
scheme, a no-rainout model would be one in which 
the fraction of condensed material follows the gas 
pressure; i.e, q t equals a constant throughout the 
cloud. In this case, the ratio of condensate to con- 
densable vapor, q c /q v , will increase with height 
above the base of the cloud. 

With the above considerations, it is straight- 
forward to calculate q c as simply the difference 
between the total number mixing ratio of condens- 
able material and the number mixing ratio of con- 
densable vapor: q c = q t — P sat /P. We therefore 
obtain the following equation for q c : 

P V 

qc '-'max ' p '-'max ' qbclow ' p • v^^J 

In Equation 21, P c , 1 denotes the saturation va- 
por pressure curve — given by P sa t , as in Equation 
1 — adjusted to resemble a condensation curve by 
applying the mixing ratio in a solar composition 
gas; i.e., P sat = q bc iow ■ P c ,i- For example, the 
dotted curve labeled "Fe" in Figure 1 shows the 
P c . i for iron, not the iron saturation vapor pres- 
sure curve. The subscript notation for P Cj i simply 
refers to the "effective" iron condensation curve as 
derived from the saturation vapor pressure curve. 
The rainout is expressed in terms of P c , i for com- 
parison with the analogous prescription for rain- 



out in a heterogeneous cloud (Equation 22). The 
quantity P in Equation 21 is the total gas pressure 
in the atmosphere. As is clear from comparing the 
thick dotted Fe line in Figure 1 with one of the gas 
pressure profiles, P decreases with height above 
the cloud base much less steeply than P c , i- The 
assumption of a constant supersaturation through- 
out the cloud, which we have made in fixing the 
value of <S max , leads to a cloud with a very low 
ratio of condensed material to total condensable 
material unless <S max 3> 0.01. 

For heterogeneously condensing clouds, we re- 
fer to the condensation curves (denoted P C) i in 
Equations 21 and 22) of forsterite and gehlenite as 
shown in Figure 1. We make the assumption for 
a heterogeneously condensing cloud that the total 
mixing ratio, q t , follows the condensation curve 
above the cloud base. This assumption is analo- 
gous to the assumption used to derive Equation 21 
for homogeneously condensing clouds that the to- 
tal mixing ratio follows the saturation vapor pres- 
sure curve. However, in a heterogeneous (or chem- 
ical) cloud, once the condensate appears in chemi- 
cal equilibrium, we argue that the solid phase will 
then be strongly favored, and hence the product 
chemical will be fully condensed. Thus, for het- 
erogeneous condensates, q v = and q c = q t . The 
rainout formula for heterogeneous species is there- 
fore quite similar to Equation 21, without the re- 
duction of q c relative to q t by the supersaturation 
factor, which is not relevant for chemical clouds: 

q c = qbciow • ~p^- (22) 

In Equation 22, the quantity P Cj i represents the 
chemical equilibrium condensation curve of the 
condensing species, as shown for forsterite and 
gehlenite in Figure 1. It is apparent from Figure 
1 that this rainout prescription sequesters most of 
the condensate within a scale height of the cloud 
base, since the condensation pressure drops off so 
steeply relative to the gas pressure. For exam- 
ple, consider the forsterite cloud shown in Fig- 
ure 1 for the brown dwarf model at T e ff = 900 K 
and g S urf = 3 x 10 5 cms -2 . At the forsterite 
cloud base, which is at about T « 2000 K and 
P « 125 bars, the mixing ratio of forsterite is just 
determined by the abundance of magnesium in the 
mixture, or q c = qbciow = 3.2 x 10~ 5 . Higher up in 
the atmosphere (in altitude), at T « 1700 K and 
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P w 70 bars, Equation 22 shows that q c , which fol- 
lows the condensation curve, is reduced by a factor 
^100 from its value at the cloud base. However, 
the gas pressure is down by only about 55%, or 
scale height. 

Equations 21 and 22 show that, in the present 
model, virtually all of the atoms in the solar- 
composition gas that are initially available for con- 
densation become sequestered in the condensate 
within a scale height of the atmosphere. For exam- 
ple, because the limiting atomic species for the for- 
mation of forsterite is magnesium, and forsterite 
is favored chemically below about 2000 K, magne- 
sium will have rained out near the forsterite cloud 
level, and no magnesium will be available for fur- 
ther condensation into other chemicals high above 
the forsterite cloud level. Rainout thus depletes 
the atmosphere of the most refractory elements at 
progressively lower temperatures and pressures. 

Each timescale characterizes the average (or 
characteristic) effectiveness of each of the compet- 
ing microphysical processes in the system. Our 
computed particle sizes thus represent the most 
probable particle size in the distribution of parti- 
cle sizes appearing at each pressure level. To the 
order of the accuracy of the particle sizes them- 
selves, therefore, the resulting particle sizes in each 
pressure level will equal the mode of the particle 
size distribution. 

We make no attempt in the present model to 
calculate the particle size distribution for each 
cloud layer. Rather, for the purposes of calculat- 
ing the opacity of clouds and using these opacities 
in SMO spectral models, we assume a functional 
form for the particle size distribution that is con- 
sistent with measurements of particle size distribu- 
tions attained in Earth's water clouds (Deirmend- 
jian 1964). We use the form of the Deirmendjian 
(1964) exponentially decaying power law employed 
by Sudarsky et al. (2000): 



2.4. Free Parameters 
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where n(r)dr is the number of particles per cubic 
centimeter having radii r — > r + dr. We employ 
the cloud model described above to determine the 
mode of the distribution, ro. 



maximum su- 
persaturation, the free parameters of the cloud 
model include the elemental abundances and the 
sticking coefficients for coagulation and coales- 
cence. The vapor mixing ratios have been cal- 
culated assuming solar (Anders & Grevesse 1989) 
abundances of the elements in the atmosphere. 

Coagulation and coalescence are not expected 
to be 100% efficient in all cases because some 
molecules stick together when they collide more 
easily than others. Therefore, we have divided the 
timescales for coagulation and coalescence, given 
in Rossow (1978) as r coag and r coa i, by efficiency 
parameters e CO ag and e coa i. 

Coalescence between solid particles is extremely 
inefficient (Rossow 1978). Coalescence will there- 
fore not be important for forsterite and gehlcn- 
ite because the solid phase is strongly favored 
once the temperature decreases sufficiently for the 
chemical to appear (see Figure 1). Therefore, for 
the hctcrogcncously condensing materials, we have 
ignored coalescence (e coa i = 0). For the species 
condensing into liquid particles (iron and water), 
we have taken e coag = 10 _1 and e coa i = 10~ 3 
(Lunine et al. 1989). We have also used e coag 
= 10 _1 for the coagulation efficiency of heteroge- 
neously condensing species. Subsequently, we will 
argue that this choice is relatively unimportant 
because coagulation and coalescence generally op- 
erate more slowly, for the atmospheres we are con- 
sidering here, than heterogeneous nucleation and 
convection. We discuss in Section 3.4 how varying 
these parameters affects the general features of the 
cloud. 

3. Model Results 

3.1. Modal Particle Sizes 

We compute particle sizes for each of four 
species on a set of atmospheric models span- 
ning the ranges T e g = 600 - 1600 K and g sur f = 
1.78 x 10 3 - 3 x 10 5 cm s~ 2 . We chose the 
most abundant condensable species in the solar- 
composition mixture: iron, forsterite, and water. 
Iron and forsterite are the most abundant high 
temperature condensates (Lunine et al. 1989). 
Although we have included one of the calcium- 
aluminum silicates, gehlenite, and these refractory 
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species do condense into clouds, their abundance 
in a solar-composition gas is lower than that of 
iron or forsterite by a factor of about ten (Lunine 
et al. 1989). Thus, the contribution of the cal- 
cium and aluminum silicates to the total opacity 
is expected to be relatively minor. Nevertheless, 
the condensation of the calcium and aluminum 
silicates is important for sequestering calcium and 
aluminum at depth beneath the photospheres of 
cooler SMOs. Moreover, the most refractory cal- 
cium and aluminum silicates can potentially serve 
as nucleation sites for the condensation of iron 
vapor via heterogeneous nucleation at low super- 
saturations. 

3.2. Trends With Brown Dwarf Gravity 
and Effective Temperature 

Figures 2,3, and 4 show the particle sizes com- 
puted by our model for a range of brown dwarf 
atmospheric profiles. Each circled point repre- 
sents, for the atmospheric profile to which it ap- 
plies, the grain size computed by our model at 
the initial cloud level (i.e., the cloud base). As 
we discuss further in Section 3.3, these are in gen- 
eral the largest modal particle sizes for the cloud, 
since particle sizes decrease with height above the 
cloud base. The model atmosphere corresponding 
to each circled point is indicated by its position on 
the diagram: gravity increases from left to right, 
whereas effective temperature increases from bot- 
tom to top in the field. These graphs say noth- 
ing about the distribution of material above the 
cloud base. They apply only to the initial con- 
densation level where the mass density of material 
is greatest. Figure 2 shows the results for iron 
grains, which condense directly from iron vapor. 
Figure 3 shows the analogous results for forsterite 
grains, which appear as a result of chemistry (i.e., 
heterogeneous condensation). Figure 4 shows the 
same for gehlenite, which also condenses heteroge- 
neously. 

Figures 2,3, and 4 show that cloud grains in hot 
brown dwarfs, for a fixed surface gravity, are sys- 
tematically larger than in cold brown dwarfs. Sim- 
ilarly, cloud particles in objects with high surface 
gravities, for a fixed effective temperature (e.g., 
older brown dwarfs), are systematically smaller 
than in objects having low surface gravities. 

The trend with surface gravity is simple to ex- 
plain: particles settle more quickly in high-gravity 
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Fig. 2. — The grain sizes of iron (Fe) at the 
cloud base for a range of brown dwarf atmospheric 
temperature-pressure profiles. The figure displays 
the modal particle sizes obtained by applying the 
model described in Section 2 to each model brown 
dwarf atmosphere. The profiles span a range of 
effective temperatures and surface gravities from 
T cff = 600 - 1600 K and g surf = 1.78 x 10 3 - 
3 x 10 5 cms -2 . Each circled point corresponds to 
a different brown dwarf model atmosphere whose 
effective temperature and gravity depend on the 
point's location in the field. The figure shows that 
high-gravity brown dwarfs will have characteristi- 
cally smaller particle sizes than low-gravity brown 
dwarfs for a given effective temperature. Simi- 
larly, hotter brown dwarfs will exhibit larger par- 
ticle sizes, for a given surface gravity, than cooler 
brown dwarfs. 
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Fig. 3. — Companion to Figure 2: shows simi- 
lar trends for forsterite (Mg 2 Si04) particle radii. 
The forsterite particles become extremely large 
(~350 /xm) in very low gravity objects. Because 
forsterite forms in abundance in the atmosphere, 
the optical effects of forsterite clouds are poten- 
tially very important, depending on the particle 
sizes. 
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Fig. 4. — The particle size trends for the calcium- 
aluminum silicate gehlenite (Ca 2 Al 2 Si07). This 
species is less abundant than forsterite by about 
a factor of ten, but it has the important effect of 
sequestering Ca and Al in condensed form below 
the photosphere of the cooler objects. 



environments. This appears to be the overrid- 
ing effect governing particle size. Although both 
the atmospheric pressure and temperature at the 
cloud base are higher in high-gravity objects, and 
the convection therefore more rapid, gravity dom- 
inates and the resulting particles are smaller. For 
the lowest gravity objects, the forsterite particles 
grow quite large. They therefore contribute only 
negligibly to the total optical depth (Sudarsky 
et al. 2000). 

The trend with effective temperature results 
from the difference in convective updraft veloc- 
ity between colder and hotter objects. Objects 
with higher effective temperatures must transport 
a higher flux of thermal energy to the surface via 
interior turbulent convection. Therefore, the up- 
drafts holding cloud particles aloft will be more 
rapid in the higher T c ff objects, resulting in larger 
particle sizes. 

3.3. Cloud Decks in Convective vs. Radia- 
tive Regions 

Brown dwarfs are almost fully convective, with 
only a thin radiative atmosphere (Basri 2000; Bur- 
rows et al. 2001). Convective uplifting sustains the 
particles against gravitational fallout. The effec- 
tiveness of this process depends on the velocity 
scale of the convection. Hence, the gradient in 
the convective velocity plays a role in determining 
particle size. Deep in the atmosphere, the brown 
dwarf is fully convective, with a relatively high 
updraft velocity (approaching 5 x 10 3 cm/s). At 
higher altitudes, the updraft speed is small, but 
not zero. Therefore, clouds at depth are convec- 
tive; clouds at altitude are quiescent. By quies- 
cent, we mean clouds formed in convectively sta- 
ble layers where the true lapse rate is shallower 
than the adiabatic lapse rate. 

The updraft speeds observed in stably strati- 
fied layers in Earth's atmosphere typically range 
from 1 — 10 cm/s. Clouds analogous to Earth's 
stratus clouds are possible in quiescent layers of 
planetary atmospheres, even though the updraft 
speeds are smaller by a factor of a hundred or 
more than the typical speeds of convective up- 
drafts (Rogers & Yau 1989). In the case of these 
stratiform clouds, small particle sizes (generally 
less than 10 /zm) lead to relatively slow sedimen- 
tation rates. Large-scale uplifting of a stable layer 
can thus replenish the cloud material lost by evap- 
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orating particles on a timescale of hours or days, 
which is rapid enough to sustain the cloud. Slow, 
large-scale updrafts, which we assume are occur- 
ring in these atmospheres, are too slow to affect 
the particles directly. Therefore, we do not include 
them as a growth process in the cloud model, as 
we have done for convective updrafts, but their 
existence is an important feature of the formation 
of quiescent clouds by maintaining the supersatu- 
rated environment necessary for particle growth. 

Particles in quiescent clouds are generally 
smaller than particles in convective clouds. Figure 
5 compares the particle sizes of condensed iron for 
a brown dwarf model having T c g- = 1500 K and 
gsurf= 5.62 x 10 4 cm s~ 2 . The modal particle 
sizes in the cloud are shown for iron under two 
contrasting assumptions: (1) the cloud is convec- 
tive (i.e., r conv = H/v conv ), and (2) the cloud is 
radiative (i.e., T con v = oo). Our calculations sug- 
gest that for this atmosphere, a convective iron 
cloud is more realistic. We thus predict the solid 
curve for the actual particle size distribution in 
the cloud deck. Figure 5 compares the particle 
sizes obtained by activating and deactivating the 
convective uplifting mechanism. 

Figure 5 shows, for this particular brown dwarf 
atmosphere, that cloud particles growing in a re- 
gion of strong convective upwelling are largest near 
the cloud base but then diminish rapidly with 
height. This decrease in size is due to the proxim- 
ity of the iron cloud base, for this particular model, 
to the radiative-convective boundary of the atmo- 
sphere, where the gradient in the upwelling ve- 
locity is steep. For convective clouds deep within 
the convective zone, where the upwelling veloc- 
ity is nearly constant, the particle sizes will be 
more constant throughout the cloud deck. How- 
ever, deep clouds will not strongly influence the 
object's emergent spectrum. 

If the cloud is radiative, the particles will not 
grow as large as they would under the influence 
of convection. Furthermore, their sizes will not 
decrease by more than a factor of two from base 
to cloud top. This result appears quite general for 
the radiative clouds we have investigated. 

The structure of the cloud decks (e.g., as shown 
in Figure 5) of brown dwarfs is crucial for com- 
puting the optical depth of the cloud. For clouds 
in which the overall particle sizes decrease sig- 
nificantly with altitude above the cloud base, it 
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Fig. 5. — The structure of an iron cloud deck for 
a particular brown dwarf at T e ff = 1500 K and 
Ssurf = 5-62 x 10 4 cm s~ 2 , assuming a constant 
1% maximum supcrsaturation during the entire 
particle growth phase. The dashed curve assumes 
that convection is not effective at sustaining parti- 
cle growth. Therefore, the particles do not grow as 
large. For this radiative (or quiescent) cloud, the 
particle sizes vary by about a factor of two within 
an atmospheric scale height. On the other hand, 
for the convective cloud, in which we allow up- 
drafts to sustain particle growth against gravita- 
tional sedimentation, the particles grow relatively 
large. Because the cloud base is slightly above 
the radiative-convective boundary of the profile, 
the updraft velocity decreases rapidly with height. 
Thus, for clouds forming near unity optical depth, 
we expect variable particle sizes from the cloud 
base to the cloud top. The cloud particle sizes 
for the convective cloud merges with the sizes ob- 
tained for the quiescent cloud once the updraft 
velocity becomes negligible. The exact cutoff for 
the cloud is somewhat arbitrary; we present a full 
scale height of both cloud decks. As we explain in 
our discussion of rainout in Section 2.3, the cloud 
is extremely thin this high above the cloud base, 
and the mole fraction of total condensing mate- 
rial, given by q t , is depleted from its value at the 
cloud base, where q t = qbciow, by many orders of 
magnitude (see Equation 21). 
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may be that the bulk of the optical depth in the 
near-infrared will be contributed not by the low- 
est cloud layers but by intermediate layers having 
significant abundances of smaller particles. 

3.4. Varying the Free Parameters 

Among the model free parameters discussed 
above, the level of maximum supersaturation is 
the most problematic. The value of iS max is par- 
ticularly difficult to calculate independently, and 
because the condensation timescale depends on it, 
<S max has a significant effect on the resulting par- 
ticle size. This effect will be manifest in quies- 
cent clouds but will be less important in convec- 
tive clouds. This is because in rapidly convecting 
regions, convective upwelling can directly sustain 
particle growth. 

Figure 6 shows, for one brown dwarf model, the 
effect of varying S max on the modal particle size of 
a water cloud. For low levels of supersaturation, 
the particles grow to a radius of ~50 ^m. In this 
case, abundant condensation nuclei raise the rate 
of heterogeneous nucleation high above the rate 
of homogeneous nucleation, and condensation can 
begin as soon as the atmosphere becomes slightly 
supersaturated. For high levels of supersaturation 
(>10%), they grow much larger, to ~150^m. In 
this clean atmosphere case, heterogeneous nucle- 
ation does not occur, and the partial pressure of 
water in the atmosphere must greatly exceed the 
saturation vapor pressure in order for condensa- 
tion to begin. We have assumed an intermediate 
value of S max = 10~ 2 for iron. This value is our 
best estimate for the <S max of vaporous iron based 
on a knowledge of the maximum supersaturations 
attained in terrestrial water clouds (Rossow 1978). 

The other free parameters of the cloud model 
are less difficult. The mixing ratios depend on the 
metallicity, which we have taken to be solar. The 
chemical equilibrium code requires the metallicity 
to be specified. Once the assumption of elemental 
composition is made, however, we can calculate 
the partial pressures of all the gases in the mixture. 
We defer a discussion of the effect of brown dwarf 
metallicity on cloud structure to future work. 

The efficiency parameters for coagulation and 
coalescence, e coag and e coa i, are potentially impor- 
tant. However, our model results show that the 
timescales for these processes, whatever their effi- 




Fig. 6. — The modal particle sizes, as a func- 
tion of iS max , for a water cloud in the upper at- 
mosphere of a cool brown dwarf: T e g = 350 K 
and g surf = 10 4 cm s~ 2 . Convection is not occur- 
ring in these outer layers. The figure demonstrates 
the difficulty in estimating particle sizes for radia- 
tive (i.e., quiescent) clouds because of the strong 
dependence of particle size on the assumed maxi- 
mum supersaturation. Water cloud particles will 
grow extremely large in a clean atmosphere; i.e., 
one in which the paucity of condensation nuclei re- 
quires high levels of supersaturation of the vapor 
in order for nucleation to begin. In this case, nu- 
cleation proceeds rapidly and the particles grow 
very large. For convective clouds, on the other 
hand, the value assumed for the maximum super- 
saturation is not so important, as the cloud parti- 
cles will grow to nearly the same size for any value 
of <S max < 100%. 
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ciency, are consistently longer than the timescale 
of growth by heterogeneous nucleation at 1% su- 
persaturation. The model suggests that coales- 
cence is much faster than coagulation for the two 
species — iron and water — in which coalescence op- 
erates efficiently. The coalescence timescale, for 
the largest particles we have grown (^500 /Um), 
can approach the timescale of heterogeneous nu- 
cleation. However, for the clouds in most of 
the atmospheres we have treated, coalescence is 
still slower than heterogeneous nucleation by a 
factor of ten or more. Indeed, unless the effi- 
ciency of coalescence is greater than one, requiring 
charged aerosols (Rossow 1978), coalescence will 
never dominate particle growth. 

4. Comparison With Other Cloud Models 
4.1. Ackerman & Marley (2001) 

Ackcrman & Marley (2001) also employ a one- 
dimensional cloud model in which cloud layers are 
treated as horizontally uniform; i.e., in a globally 
averaged sense. The Ackerman & Marley (2001) 
cloud model arrives at particle sizes by balanc- 
ing the upward transport of condensable mate- 
rial due to convective updrafts with the down- 
ward sedimentation of condensate. In the case in 
which the convective timescale, r conv , dominates 
over the other growth timescales in the present 
cloud model, which will occur in vigorously con- 
verting regions of the atmosphere, the results of 
the present cloud model are expected to produce 
very similar results to the Ackerman & Marley 
(2001) treatment. This is because our formula- 
tion then reduces to a similar comparison between 
the characteristic velocities of turbulent convective 
updrafts and the terminal velocity of particles. 

In their treatment of rainout, Ackcrman & 
Marley (2001) introduce an additional parameter, 
/rain, which characterizes the vertical distribution 
of condensate in the cloud. The / ra in parameter, 
as Ackerman & Marley (2001) state, is difficult 
to calculate from basic principles. They leave it 
as an adjustable parameter in their model. It de- 
pends on the mass-weighted sedimentation veloc- 
ity of the cloud droplets and the convective veloc- 
ity scale. The difficulty in calculating / ra j n lies in 
the complexity of modeling fully the nature of the 
convection within the cloud. This problem has yet 
to be successfully attacked for brown dwarf atmo- 



spheres. 

If it can be computed for these convective 
clouds, however, knowledge of / ra i n can potentially 
provide a more realistic measure of the height and 
distribution of material in the cloud. For quies- 
cent clouds, / ra i n does not play a role. Acker- 
man & Marley (2001) have assumed / ra i n = 2-3 
in their calculations. The cloud decks they derive 
are somewhat more compact than those produced 
in the present cloud model, although both pre- 
scriptions lead to rapid depletion of cloud material 
within a fraction of an atmospheric scale height 
above the condensation level. 

The particle sizes we obtain compare favorably 
with the sizes predicted by the Ackerman & Mar- 
ley (2001) model. They predict modal particle 
radii in the intermediate size range of 40-80 ^m 
for both iron and silicate grains, in good agree- 
ment with the particle size ranges shown in Fig- 
ures 2, 3, and 4. 

4.2. Helling et al. (2001) 

Helling ct al. (2001) study the onset of cloud 
particle growth via acoustic waves. They show 
that small (1 — 10 fim in radius) sized particles 
can nucleate rapidly, normally within a few sec- 
onds. These values are consistent with our par- 
ticle sizes in radiative regions under the assump- 
tion of very low supersaturations (5 max <C 10~ 2 ). 
Our larger particles form from sustained particle 
growth by convective uplifting. The grains gener- 
ated in (Helling et al. 2001) are not grown to the 
maximum size allowed gravitationally, which we 
assume can occur. 

The predictions of the two models in terms of 
particle radii will be in qualitative agreement in 
the absence of the r conv timescale employed in 
our model. Nevertheless, we emphasize caution in 
comparing directly the results of these two models. 
They employ vastly different physical approaches 
and are therefore not expected to agree in many 
cases, even qualitatively. The model of Helling 
et al. (2001) attempts a detailed simulation of 
particle growth in brown dwarf atmospheres, in- 
cluding the complicated effects of hydrodynamics. 
The goal of the present cloud model, rather, is 
not to directly simulate the dynamics of particle 
growth but to develop a computationally econom- 
ical cloud model that can be incorporated easily 
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into spectral synthesis models. 



Discussion 



5.1. Cloud Opacity 

The effects of clouds on the emergent spectra of 
substellar objects depend strongly on grain sizes. 
We derive the wavelength-dependent absorption 
and scattering opacities of grains with a full Mie 
theory approach, utilizing the formalism of van de 
Hulst (1957) and Dcirmendjian (1969). 

Figure 7 shows the results of such a calcula- 
tion for forsterite grains in a brown dwarf at- 
mosphere (T cff = 1500 K, g surf = 10 5 cm s~ 2 ) 
based on our cloud model results and optical con- 
stants from Scott & Duley (1996). Also shown 
are the results for a larger grain size distribution 
peaked at 50 /um, as well as for a size distribu- 
tion peaked at 0.1 yum, which is representative of 
an interstellar particle size (Mathis et al. 1977) 
often assumed to be appropriate for substellar ob- 
jects (Allard et al. 2001; Barman et al. 2001). 
For comparison, we have plotted the atomic and 
molecular gaseous opacities (Burrows et al. 2001, 
and references therein) within the cloud region. 
The substantial absorption and scattering differ- 
ences between the size distributions, and relative 
to the gaseous absorption, convey the importance 
of proper cloud modeling. 

5.2. Effect on Spectra 

Figure 8 depicts the effects of our modeled 
forsterite cloud on the emergent spectrum of a 
brown dwarf (T e ff = 1500 K, g surf = 10 5 cm s~ 2 ). 
This model atmosphere was obtained using the 
self -consistent stellar atmosphere code, TLUSTY 
(Hubeny, I. 1988; Hubeny & Lanz 1995). The base 
of the forsterite cloud resides at approximately 4 
bars and 1800 K, which is the highest temperature 
in this atmosphere at which forsterite grains can 
form. 

We made two simplifying approximations in 
producing the cloudy model spectrum shown in 
Figure 8. First, a cutoff for the cloud deck was de- 
sirable to facilitate the calulation of optical depth. 
We chose this cutoff to be one scale height above 
the cloud base, which is the atmospheric pressure 
e-folding distance. Since the remaining cloud- 
forming material is highly depleted of condensate a 



cloud model result (11-micron grains): 
absorp. scat. 

50-micron grain distribution: 
absorp. scat 



0.1 -micron grain distribution 
absorp. scat. 




wavelength (microns) 

Fig. 7. — Scattering and absorption opacities of 
forsterite are compared with gaseous atomic and 
molecular absorption at (1600 K, 2.5 bar). Our 
modeled modal particle size, 11 /*m, in a brown 
dwarf (T cff = 1500 K, g surf = 10 5 cm s~ 2 ) is con- 
trasted with both a larger size distribution peaked 
at 50 /im and a smaller size distribution peaked at 
0.1 /jm (representative of an interstellar grain size 
assumed by some researchers). In all cases, a func- 
tional form of the distribution about the modal 
size is used (Dcirmendjian 1964). 
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Fig. 8. — The emergent spectrum of a cloudy 
brown dwarf (T e ff = 1500 K, g surf = 10 5 cm s~ 2 ) 
is compared with a cloud-free model of the same 
effective temperature and surface gravity. In the 
cloudy model, the base of the forsterite cloud re- 
sides at approximately 4 bars (1800 K). For the 
synthetic cloudy spectrum shown, we took the 
cloud deck to be a scale height thick. This choice 
will not affect the shape of the model spectrum, 
however, so long as the cloud deck is more than 
about a third of a scale height thick, since the 
cloud is very tenuous and thus optically thin high 
above the cloud base. The particle size distribu- 
tion was chosen to be the same throughout the 
cloud deck for simplicity. We used a Deirmendjian 
(1964) distribution centered around 11 /im, which 
was the typical particle size in the forsterite cloud 
deck as computed by the cloud model. Within the 
cloud region, the full solar abundance of magne- 
sium has been condensed into forsterite grains. 



scale height above the condensation level (see Sec- 
tion 2.3, Equation 22), the cloud in this region is 
optically thin. This simplification is therefore not 
a concern, and the spectral model shown in Figure 
8 incorporates the full opacity of the cloud. 

Second, owing to difficulties in computing the 
Mie theory scattering and absorption opacities it- 
eratively, we chose to simplify the radiative trans- 
fer problem by using a uniform modal particle ra- 
dius of 11 ^m throughout the cloudy region. That 
is, at every atmospheric pressure level for which 
we computed the opacity of particles, we employed 
the Deirmendjian (1964) particle size distribution, 
given by Equation 23, with value of r equal to 
1 1 ^m. The scattering and absorption opacities of 
a cloud of 1 1 /im forsterite particles are shown in 
Figure 7. This value represents the typical particle 
size in the cloud deck, accounting for the decrease 
in density of cloud material with height above the 
cloud base (i.e., a number density weighted av- 
erage of the particle sizes predicted by the cloud 
model). 

A cloud-free model of the same effective tem- 
perature and gravity is plotted for comparison. 
Within the B (~ 0.45 fim) and Z (<~ 1 fim) bands, 
the emergent flux is lowered by as much as one 
dcx. However, the strong absorption by the wings 
of the sodium and potassium resonance lines is 
mitigated somewhat due to the clouds. Also of 
interest are the differences in the J (~ 1.25 /im), 
H (~ 1.6 ^m), and K (~ 2.2 ^m) infrared bands. 
The presence of clouds reduces the emergent flux 
in the otherwise relatively clear J and H bands, 
allowing more flux to escape between these bands 
and at longer wavelengths. In fact, the J — K 
colors differ by over 1.5 magnitudes between the 
cloudy and cloud-free models. 

The cloudy spectrum presented in Figure 8 is 
intended to demonstrate the potential importance 
of clouds as an opacity source in SMO atmo- 
spheres. We defer to future work the problem 
of more realistically incorporating the opacity of 
clouds into spectral synthesis models. 

5.3. Coupling Clouds With SMO Atmo- 
sphere Models 



Developing fully self-consistent atmosphere 
models with clouds is potentially problematic be- 
cause clouds perturb the radiative balance of the 
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upper atmosphere. Scattering and absorption of 
radiation by the cloud causes the temperature- 
pressure structure of a cloudy atmosphere to devi- 
ate from the structure of a cloud-free atmosphere 
having the same effective temperature and sur- 
face gravity. The extent of the effect is strongly 
dependent on the location of the cloud in the 
atmosphere and the vertical variation of cloud 
particle sizes and number densities. Calculating 
an atmospheric profile that is in radiative equi- 
librium including both the gaseous opacity and 
the opacity of the cloud is straightforward. The 
potential inconsistency arises from the fact that 
the nature of the cloud itself, including the cloud 
base location and the particle sizes, depends in 
general on the temperature-pressure profile of the 
atmosphere, as we show in Section 2. This is a 
problem faced by all research groups attempting 
to include clouds completely self-consistently into 
SMO spectral synthesis models. 

We explored this problem quantitatively by 
comparing the temperature-pressure profile of two 
atmosphere models, a "cloudy" model atmosphere 
and a "cloud-free" model atmosphere, at the same 
effective temperature and surface gravity. The 
cloud-free model used for this test was obtained 
using TLUSTY, our self-consistent model atmo- 
sphere code (Hubeny, I. 1988; Hubeny & Lanz 
1995), by including only the opacity from gaseous 
atomic and molecular absorption. We then ap- 
plied the cloud model described in Section 2 to this 
cloud-free atmosphere to obtain the distribution 
of particle sizes and densities of the forsterite cloud 
predicted to form near the object's visible surface. 
For simplicity, we took a sensible average of the 
particle sizes output by the cloud model and cal- 
culated, for a cloud of uniform modal particle size, 
the absorption and scattering of radiation by the 
particles using the Mie theory approach outlined 
in Section 5.1. We then incorporated this opac- 
ity back into the atmosphere code to recompute 
the temperature-pressure structure, thus obtain- 
ing the cloudy model atmosphere. We found the 
cloudy atmosphere to be hotter by several hun- 
dred degrees than the cloud-free atmosphere at 
the same pressure, a significant change. This is 
the back-warming effect of the cloud alluded to 
previously. 

We then applied the cloud model to the new 
cloudy atmosphere to see whether the forsterite 



grain sizes varied significantly as a result of the 
change in the temperature-pressure profile. We 
found that the particles did change in size: their 
radius increased by about a factor of three. That 
is, the forsterite cloud looks quite different for the 
perturbed atmosphere model than it does for the 
original cloud-free model. This change is a re- 
sult of the fact that the forsterite cloud straddles 
the convective-radiative boundary of the atmo- 
sphere. Thus, in the original model, the forsterite 
cloud appears in a convectively stable region of 
the brown dwarf, but in the perturbed model, the 
forsterite cloud forms in the convective zone of the 
atmosphere. Thus, the particles in the perturbed 
model came out larger than the particles in the 
original model (Figure 5 shows how particles in 
radiative regions are systematically smaller than 
particles forming in convective regions). 

We performed a similar test for iron clouds. Un- 
like for the forsterite cloud, we found only a small 
change in the particle sizes — a decrease of about 
15% in radius — resulting from the perturbation in 
the atmospheric temperature-pressure structure. 
The iron cloud forms deeply enough in the con- 
vective region of the brown dwarf studied to not 
be strongly affected by the heating due to the 
cloud. The degree of inconsistency between cloud- 
free and cloudy profiles therefore depends on the 
details of the calculation itself. In many cases, 
the discrepancies will not be a major concern, but 
they can be particularly large if the condensation 
curve happens to intersect the atmospheric profile 
near the convective-radiative boundary. 

We emphasize that both the "cloud-free" and 
"cloudy" model atmospheres described above sat- 
isfy the radiative equilibrium boundary condition, 
and in that sense, the profiles are self-consistent. 
The remaining inconsistency arises from not ac- 
counting properly for possible adjustments to the 
opacity of the cloud when the particle sizes change 
after the profile is perturbed. A future chal- 
lenge will be developing the machinery to generate 
model atmospheres that incorporate clouds fully 
self-consistently. 

6. Conclusions 

6.1. General Results 

We have addressed the condensation and subse- 
quent growth of cloud particles in the atmospheres 
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of brown dwarfs. We present optimal particle sizes 
for three abundant species — iron, forsterite, and 
gehlenite — for a broad range of brown dwarf effec- 
tive temperatures and surface gravities. 

High-gravity brown dwarfs exhibit clouds with 
typical particle sizes in the 5 — 20 ^m range. 
The particles grow much larger, however, in low- 
gravity objects, often greater than 100 ^m in ra- 
dius. We discovered a similar trend with effective 
temperature: hot brown dwarfs have characteris- 
tically larger particle sizes than cool brown dwarfs 
because of the increased energy flux that must be 
transported by convection. The distribution of 
cloud particle sizes depends strongly on the atmo- 
spheric parameters, and it is therefore unrealistic 
in spectral models to assume a single particle size 
distribution for the entire class of SMOs. 

We demonstrate that particle size crucially af- 
fects the optical depth of the cloud. Unlike clouds 
having a particle size distribution centered at 
0.1 /urn, these cloud decks do not dominate the 
opacity. Rather, they smooth the emergent spec- 
trum and partially redistribute the radiative en- 
ergy (see Figure 8). 

6.2. Improved Cloud Models: Morphol- 
ogy and Patchiness 

A complete theory of brown dwarf and giant 
planet atmospheres will require detailed modeling 
of clouds. Although obtaining plausible particle 
sizes is a good start, we have been unable to say 
anything yet about the morphology or patchiness 
of clouds. For example, on Jupiter, clouds appear 
in banded structures that vary latitudinally, while 
Neptune and Uranus show less latitudinal banding 
in their surface cloud patterns. The structure and 
patchiness of clouds on brown dwarfs is not known, 
though the putative variability reported by Bailer- 
Jones & Mundt (2001) suggests that cloud patches 
are not necessarily static. 

A simple approach to understanding surface 
cloud distributions is the moist entraining plume 
model. Such a model has been put forth by Luninc 
ct al. (1989) based on the work of Stoker (1986) 
for Jupiter's equatorial plumes. The model gen- 
erates plumes from atmospheric material whose 
buoyancy is increased by the release of latent heat 
from condensation. 

More elaborate calculations would involve 



three-dimensional modeling of the general atmo- 
spheric circulation. Such an approach would be 
quite difficult to implement, especially without ac- 
curate meteorological data for SMO atmospheres. 
General circulation models would also be too in- 
tensive computationally to directly incorporate 
into SMO spectral models. 

6.3. Future Work in Spectral Synthesis 

The results of the current work will be use- 
ful in developing more elaborate spectral models 
of substellar atmospheres. We plan to follow up 
the present work with an exploration of the spec- 
tral effects of a variety of cloud compositions and 
distributions. Although Figure 8 shows the gen- 
eral effects of introducing a forsterite cloud on the 
spectrum of a brown dwarf, we have not yet var- 
ied the particle size distribution from the base to 
the top of the cloud. A large drop in the particle 
sizes with height above the cloud base will pro- 
duce greater optical depth than a cloud deck hav- 
ing uniformly large particles. It is also likely that 
clouds of more than one species will form near the 
photosphere simultaneously A major future chal- 
lenge of this field will be to model the time de- 
pendent dust formation processes of all the major 
condensable species together, including the pos- 
sible formation of core-mantle grains, which are 
composed of multiple chemical phases, and then to 
incorporate them fully self-consistently into sub- 
stellar model atmospheres. 

We thank Ivan Hubeny, Chris Sharp, Jonathan 
Fortney, Jason Barnes, Adam Showman, Bill Hub- 
bard, Jim Liebcrt, and the referee for insightful 
discussions and advice. This work was supported 
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